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Abstract 

The B20 compound MnSi exhibits a skyrmion phase in a certain range of temperature (T) and 
external magnetic field (H) in addition to the helimagnetic ordering. In this work, the magnetic 
entropy change (ASj,) in the single crystal helimagnet MnSi has been investigated by the scaling 
analysis. The ASm(T, H) curves in higher field range (2 kOe < H < 5.6 kOe) collapse onto an 
universal curve after scaling, which is independent of T and H. Based on this universal curve, the 
ASm(T, H) (excluding the skyrmion entropy change) at any T and H can be calculated. On the 
other hand, a small valley is found on the ASm(T) curve at ~ 28.5 K in the field range from ~ 1 kOe 
to ~ 2 kOe, which is resulted from the skyrmion phase. Subsequently, the magnetic entropy change 
of the skyrmion phase [AS*** (T, H)] can be derived by subtraction of the calculated AS¢4!(T, H) 
from the experimental AS%7” (T, H). The obtained results show that AS SEX reaches the maximum 


at ~ 1.6 kOe, which indicates the highest number of skyrmion vortices at that field. 
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I. INTRODUCTION 


Recently, the skyrmion ordering has triggered great interest because of its great potential 
application in spintronics storage device due to the topological property [1 2], nanometric 
size [8], and current-driven motion |4\|5,'6). The skyrmion phase is a topologically stable spin 
vortex-like texture existing in a certain range of temperature and magnetic field |7|. The 
compounds exhibiting a skyrmion phase, such as FeGe [8], MnSi [9], Feı_,Co,Si Li], 
CusOSeO3 (13), etc, have been extensively and intensively investigated. Especially, 
much interest has been drawn on the B20 compound MnSi due to the additional weak 
itinerant ferromagnetism [14], non-Fermi-liquid behavior [15, 16], and magnetic quantum 
phase transition |17|. The cell of MnSi has cubic symmetry with space group P213. The lack 
of centrosymmetry results in a weak Dzyaloshinskii-Moriya (DM) spin-orbit coupling. The 
DM spin-orbit coupling competes with the much stronger ferromagnetic exchange, which 
finally causes a long modulation period of a helimagnetic ground state [18]. It has 
been demonstrated that MnSi exhibits a first-order phase transition into the helimagnetic 
ordering at 30.5 K, which can be suppressed into a second one by an external field where a 
tricritical point was induced 23]. When an external magnetic field is applied 
above a threshold value, the helimagnetic state transforms into a conical phase |7|. The 
chiral skyrmion phase emerges just below Tc in a narrow magnetic field range due to a 
unique stabilization provided by the DM interactions |1]. In the bulk material of MnSi, the 
skyrmion phase emerges at ~ 28 K in the field range from 1.55 kOe to 2.3 kOe [23]. It has 
been proved that the skyrmion phase can be stabilized by reducing the dimension of the 
materials |24) |25). 

As we know, a magnetic ordering phase transition will result into a change of the magnetic 
entropy. Recently, a scaling analysis of the magnetic entropy change correlated to critical 
exponents has been proved an effective method to study the magnetic entropy change |26} 27, 
28), 29]. In this work, the magnetic entropy change (AS),;) in the single crystal helimagnet 
MnSi has been investigated by the scaling analysis. The ASy(T) curves exhibit a valley 
at 28.5 K in the magnetic range from ~ 1 kOe to ~ 2 kOe, which just corresponds to the 
formation of the skyrmion phase. The ASm(T, H) curves for higher fields collapse onto an 
universal curve after scaling, which is independent of the temperature T and the applied 


external field H. Based on this universal curve, ASy(T, H) (excluding the skyrmion entropy 


change) at any T and H can be calculated. Subsequently, the magnetic entropy change of the 
skyrmion phase [AS*"*(T, H)] can be derived by subtraction of the calculated AS“ (T, H) 
from the experimental AS% (T, H). It is found that AS*?* reaches the maximum at ~ 
1.6 kOe, which indicates the highest number of skyrmion vortices in MnSi at this magnetic 


field. 


II. EXPERIMENT 


A single crystal sample of MnSi was prepared by the Czochralski method [30]. The spec- 
imen quality and physical properties were checked and described elsewhere |22|. The mag- 
netization was measured using a Quantum Design vibrating sample magnetometer (SQUID- 
VSM). The no-overshoot mode was applied to ensure a precise magnetic field. The applied 
field was relaxed for 2 minutes before the data collection. For the measurement the sam- 
ple was processed into slender ellipsoid shape, and the magnetic field was applied along 
the longest axis to minimize the demagnetizating field. To make sure that each curve was 
initially magnetized, the isothermal magnetization measurement was performed after the 
sample was heated well above To for 10 minutes and then cooled under zero field to the 
target temperatures. The magnetic background was carefully subtracted. The applied mag- 
netic field H, has been corrected into the internal field as H = H, — NM (where M is 
the measured magnetization and N is the demagnetization factor) [31]. The corrected H 
was used for the calculation of the magnetic entropy change. All the magnetic units and 


formulas concerned are given in Gauss system of units. 


II. RESULTS AND DISCUSSION 


Figure [1] depicts the temperature dependence of magnetization [M(T)| for MnSi under 
zero-field-cooling (ZFC) and field-cooling (FC) under an applied magnetic field H = 0.1 
kOe, where the inset shows that for H = 2 kOe. The M(T) curves for 0.1 kOe indicates 
a magnetic ordering transition at Te, which corresponds to the well-known paramagnetic- 
helimagnetic (PM-HM) transition in MnSi [33]. The transition temperature To ~ 31 
K is determined from the dM/dT curve (the right axis in the inset). However, for M(T) 


for H = 2 kOe, the magnetic transition is attributed to ordering from a paramagnetic to 


conical spin state, where the helimagnetic ordering is polarized into conical ordering by the 
external magnetic field [7]. 

In order to investigate the magnetic entropy change, the initial isothermal magnetization 
M(H) around To for MnSi was measured, as shown in Fig. |2|(a). For the M(H) curve far 
below Te, it increases linearly with the increase of H below 5.6 kOe, and become almostly 
independent of H when exceeding 5.6 kOe. This result indicates the spin ordering behavior in 
this system. Below Tọ, the spins are polarized gradually when H below 5.6 kOe, and becomes 
completely polarized when field exceeding 5.6 kOe. As it is well known, the magnetic entropy 


change (AS),) can be calculated as [34]: 


H OCT) 
AS (T, H) = a |) daH 1 
sam- | (FP), (1 
With Maxwell’s relation: [3S(T, H)/OH|r=[OM(T, H)/OT]r |35|, ASm(T, H) is written as: 
HA /OM(T, H) 
ASm(T, H) = ———— |) dH (2) 
0 oT H 
Subsequently, ASy can be calculated from experimental data as [36]: 
= (Mi — Misi) x, 
|ASul = >, | nn |^ (3) 


where M; and M;.ı are the experimental data of the magnetization at T; and T;+ı under a 
magnetic field H;. Figure [2] (b) shows the ASm(T, H) around Tc, which is obtained from 
the isothermal M(H) curves. It can be seen that all curves almost exhibit a unimodal shape. 
However, there is an anomaly at ~ 28.5 K for the AS), just below Tg. The inset of Fig. 
2] (b) gives typical ASy(T) curves for H = 1.2 and 1.8 kOe, both of which exhibit small 
valleys at ~ 28.5 K. This valley on ASm(T) curve exists in the field range from ~ 1 kOe to 
~ 2 kOe, and disappears gradually when the applied field exceeds 2 kOe. As it is known, 
MnSi exhibits a skyrmion transition at ~ 28.5 K in the field range from ~ 1.55 kOe to ~ 
2.3 kOe [23]. The skyrmion transition should result into a change of magnetic entropy due 
to the ordering of the spins. Therefore, the valley on ASm(T) curve should be attributed to 
the entropy change of skyrmion transition. Actually, a peak caused by the skyrmion lattice 
was found in the specific heat measurement, which confirms the magnetic entropy change 
of the skyrmion phase observed here [23]. 

The ASu(T, H) can be scaled onto a phenomenological universal curve based on scal- 


ing relation 28|. According to scaling laws, the experimental ASm(T, H) should 
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collapse onto two independent universal curves below and above Tc with proper reference 
temperatures and normalization. The ASm(T) curves for different H were normalized as 
\ASu(T)/ASY"(T)| on the longitudinal coordinate, where the maximum of all curves are 
located on single point on the new normalized longitudinal coordinate. The horizontal co- 
ordinate is marked as the rescaled temperature 6, which is defined as [28]: 


ee 


p= 2 
Tel; 


(4) 


where T, is the reference temperature. By choosing proper T,, the equivalent points of the 
different experimental curves should collapse onto the same point on the universal curve 
28|. The experimental ASm(T, H) should collapse onto two independent universal curves 
below and above Tọ respectively. The combination of the two universal curves below and 
above To just form a whole universal magnetic entropy change curve. After scaling and 
normalizing, the universal curve would be to impose the position of T, at a value of 6 = 1. 

Theoretical investigations have demonstrated that it is not necessary to use two reference 
temperatures below and above Tc but just a single one for a single magnetic phase material. 
However, in the case of a large demagnetization factor or the existence of multiple magnetic 
phases, the use of two reference temperatures is necessary. In this case, it is noticed that 
ASm(T) for different H present asymmetry. In present experiment, the demagnetization 
energy has been carefully excluded as aforementioned. However, asymmetry of ASy(T) 
still appears, which may be caused by the multiple magnetic phase transitions and magnetic 
fluctuation 39|. Therefore, two reference temperatures should be used, which are 


defined as [28]: 
9- 0- = (To —T)/(In — Te), T < Te 6) 
6, = (T-Te)/(Ta — Te), T > Te 
where T,ı and T,2 are the reference temperatures below and above Te, respectively. Gen- 
erally, temperatures at the full-width-at-half-maximum (FWHM) are chosen as the T,ı and 
T,2, where T, corresponds to the one below To while T, to the one above Te |29}. Then, in 
the new scaled coordinate systems, 0 = +1 just corresponds to the positions of full-width- 
at-half-maximum (i.e. T,ı and T,2). 
The skyrmion phase has been absolutely suppressed by the external magnetic field of 
~ 2.3 kOe. Thus, there is no skyrmion magnetic entropy change for ASm(T, H) for a 
field exceeding 2.3 kOe. The scaling of ASy(T, H) for higher fields exceeding 2.3 kOe 
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can generate an universal curve, which is independent of T and H. The phenomenological 
construction of the universal scaling of the experimental ASy(T, H) is plotted in Fig. 3} 
The ASy is normalized as |[ASm/AS 7|, and the temperature is plotted as a rescaled 
temperature 0. One can see that all experimental data collapse onto a single universal 
curve, except a small dispersion just below To which may be due to the strong magnetic 
fluctuation approaching Te. Based on the universal scaling, an average curve is obtained, 
as shown as the red solid curve in Fig. 3, Obviously, the obtained average curve is universal 
and independent of the external measurement conditions, which is just determined by the 
intrinsic magnetization. According to the obtained universal curve, ASy(T, H) at any T 
and H can be calculated. The inset of Fig. 38 shows the comparison of the calculated and 
experimental curves for the typical ASm(T) for H = 1.2 kOe, where the calculated result 
agrees well with the experimental one except for that around 28.5 K. 

This universal curve can be extrapolated to low fields to yield a calculated AS“ (T, H). 
However, the entropy change of the skyrmion phase is not included by this calculation 
method. Therefore, through the comparison of the calculated AS), and the experimental 
curve, the skyrmion magnetic entropy change can be derived by subtracting the AS¢4!(T, H) 
from the experimental AS% (T, H). It can be seen that the calculated curve fits the exper- 
imental one very well except for the temperature range around 28.5 K. The discrepancy is 
attributed to the magnetic entropy caused by the skyrmion phase transition. Subsequently, 


the magnetic entropy change of the skyrmion (AS*"*) can be obtained as: 
AS (T)| = [ASHP (T) - ASH#(T) (6) 


Figure/4 (a) shows the obtained |AS'*"* (T)| for different H. It can be seen that |AS*?"* (T)| 


reaches the maximum at ~ 28.3 K. Figure [4 (b) plots the maximum of |AS**(T)| vs 


H, which exhibits that the maximum of the ASP 


subsequently decreases with H. The As reaches the maximum at H ~ 1.6 kOe. As 


(T)| firstly increases with H, and 


it is known, Ao originates from the magnetic entropy change of the skyrmion phase 
transition. Thus, the maximum of the AS**"* corresponds to the most drastic point of the 
skyrmion transition. Therefore, the results of the magnetic entropy change indicate that 
the maximal number of skyrmion vortices in MnSi is reached at H ~ 1.6 kOe, which is in 


agreement with the results of specific heat and detailed magnetic measurements [7, 23]. 


IV. CONCLUSION 


In summary, the magnetic entropy change AS), of MnSi has been investigated by the 
scaling analysis. The ASy(T) curves in the field range from ~ 1 kOe to ~ 2 kOe exhibit 
small valleys at ~ 28.5 K, which correspond to the magnetic change caused by the skyrmion 
phase transition. The ASm(T, H) curves in the higher field range collapse onto an universal 
curve after scaling, which is independent of T and H. Subsequently, magnetic entropy 
change of the skyrmion phase [AS*** (T, H)] can be derived by subtraction of the calculated 
AS (T, H) from the experimental AS$” (T, H). It is demonstrated that the AS" reaches 
the maximum at ~ 1.6 kOe, which indicates that the maximal number of skyrmion vortices 


in MnSi is reached at H ~1.6 kOe. 
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FIG. 1: (Color online) Temperature dependence of the magnetization [M(T)] for MnSi for H = 
0.1 kOe [the inset shows M(T) (left axis) and dM/dT (right one) for H = 2 kOe]. 
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FIG. 2: (Color online) (a) Initial magnetization M(H) around Tc for MnSi (the red curve corre- 
sponds that at 29 K); (b) magnetic entropy change (ASm) around Tc for different fields (the inset 


shows the typical curves for H = 1.2 and 1.8 kOe). 
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FIG. 3: (Color online) Scaling of ASm(T, H) curves for MnSi (the inset shows the comparison of 


the calculated and the experimental curve for H = 1.2 kOe). 
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FIG. 4: (Color online) (a) Magnetic entropy change of skyrmion transition [|AS 


kyX 
vr (T)|] for 
different fields as a function of temperature; (b) maximum of the IA S$9X (T)| vs H (the curves 
are guided on eye). 
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